Theory of spinor Fermi and Bose gases in tight atom waveguides 
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Divergence-free pseudopotentials for spatially even and odd-wave interactions in spinor Fermi 
gases in tight atom waveguides are derived. The Fermi-Bose mapping method is used to relate the 
effectively one-dimensional fermionic many-body problem to that of a spinor Bose gas. Depending on 
the relative magnitudes of the even and odd-wave interactions, the A^'-atom ground state may have 
total spin S = 0, S = N/2, and possibly also intermediate values, the case S = N/2 applying near 
a p-wave Feshbach resonance, where the A'^-fermion ground state is space-antisymmetric and spin- 
symmetric. In this case the fermionic ground state maps to the spinless bosonic Lieb-Liniger gas. 
An external magnetic field with a longitudinal gradient causes a Stern-Gerlach spatial separation of 
the corresponding trapped Fermi gas with respect to various values of Sz- 
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I. INTRODUCTION 

When an ultracold atomic vapor is placed into an atom 
waveguide with sufficiently tight transverse confinement, 
its two-body scattering properties are strongly modified. 
This occurs in a regime of low temperatures and den- 
sities where transverse oscillator modes are frozen and 
the dynamics is described by a one-dimensional (ID) 
Hamiltonian with zero-range interactions [1, 2], a regime 
which has already been reached experimentally [3-5] . In 
addition, a regime with chemical potential fj, less than 
transverse level spacing hLU± but ksT > hLu± has been 
achieved [6, 7]. We assume herein that both /i < huj± and 
ksT < hLUj_ as in [3-5]. Nevertheless, virtually excited 
transverse modes renormalize the effective ID coupling 
constant via a confinement-induced resonance, as first 
shown for bosons [1, 8] and recently for spin-polarized 
fermionic vapors by Granger and Blume [9]. 

The dynamics of an optically trapped Fermi gas is 
richer than that of a magnetically trapped one, since the 
spin is not polarized. A Fermi-Bose mapping first used 
to solve the ID hard-sphere Bose gas [10, 11] was re- 
cently shown [12] to provide an exact duality between 
effective zero-range ID fermionic and bosonic interac- 
tions and applied to spin-polarized Fermi gases [9, 13]. 
It will be shown here that this mapping and Fermi-Bose 
duality also hold for spinor Fermi gases. This mapping 
will be exploited to reduce the degenerate spatially an- 
tisymmetric fermionic ground states to that of a spin- 
less Bose gas, which, in the ID, zero-range interaction 
regime, is the Lieb-Liniger model which is exactly solu- 
ble in the absence of longitudinal trapping [14, 15] and 
well approximated by a local equilibrium approach in the 
trapped case [16]. It will be shown that in the presence 
of longitudinal trapping, this fermionic ground state can 



be Stern-Gerlach spatially decomposed by a longitudi- 
nal magnetic field gradient into components with various 
values of total longitudinal spin. 



II. TWO-BODY PROBLEM FOR FERMIONS IN 
A TIGHT WAVEGUIDE 

Consider first the three-dimensional two-body scat- 
tering problem for spin-i fermionic atoms. There are 
both s-wave scattering states, which are space symmet- 
ric and spin antisymmetric with spin eigenfunctions of 
singlet form -^(Ti — it), as well as p-wave scattering 
states which are space antisymmetric and spin symmet- 
ric with spin eigenfunctions of triplet form ft or || or 
-^(tJ. + it)- s-wave scattering cannot occur in a spin- 
polarized Fermi gas, but it is usually dominant in a spinor 
Fermi gas since p-wave spatial antisymmetry suppresses 
short-range interactions. However, both s-wave and p- 
wave interactions can be greatly enhanced by Feshbach 
resonances [17, 18]. Assume until further notice that the 
Hamiltonian does not depend on spin. Then the spin de- 
pendence of wave functions need not be indicated explic- 
itly and they can be written as the sum of spatially even 
and odd parts ipe and ■0o. When such an atomic vapor 
is confined in an atom waveguide with tight transverse 
trapping, the dynamics becomes effectively ID [1]. The 
effective ID interactions are determined by ID scatter- 
ing lengths alj^ for spatially even waves ipdz) = ipe{—z) 
related to 3D s-wave scattering and spatially odd waves 
tpoiz) = —ipoi—z) related to 3D p-wave scattering; here z 
is the relative coordinate zi — Z2 for ID scattering. These 
determine the A: ^ behavior just outside the range zq 
of the interaction: 

i^eizo) = -'/'e(-^o) = -{alo - 2o)"V-'e(±2o) 

^o{zo) = -V'o(-zo) = -(a?n - 2o)V'o(±2o) • (1) 
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is a known [1] function of the 3D s-wave scattering 
length Qs, and is a known [9, 19] function of the 3D p- 
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wave scattering volume = = — linife^o taii(5j,(fc)/fc^ 
[20]: 

Ki.)-' = ^[l-(a./a^)|C(l/2)|]-^ 

a_|_ 

alo = ^[1 + 12(V4)IC(-1/2,1)|]-^ (2) 

"_L 

where aj_ — y^h/ fiuj^ is tlic transverse oscillator length 
for the relative motion, /x is the effective mass, C(l/2) = 
—1.460 • • • is a Riemann zeta function, and ({—1/2, 1) = 
-C(3/2)/47r = -0.2079- •• is a Hurwitz zeta function 
[21]. 

In the zero-range limit Zq approaches 0+ and —zq ap- 
proaches 0—, and Eqs. (1) reduce to "contact condi- 
tions" which relate a discontinuity in at contact z = 
to Ve(0), and a discontinuity of ipo to V'o(O)- Although 
a discontinuity in the derivative is a well-known conse- 
quence of the zero-range delta function pseudopotential 
and plays a crucial role in the solution of the Lieb-Liniger 
model [14] , discontinuities of tp itself have received little 
attention, although they have been discussed previously 
by Cheon and Shigehara [12] and arc implicit in the re- 
cent work of Granger and Blume [9] . For an odd wave tpo 
the discontinuity 2'iIj{0+) is a trivial consequence of anti- 
symmetry together with the fact that a nonzero odd- wave 
scattering length cannot be obtained in the limit zq ^ 
unless i/'o(0±) ^ 0. These discontinuities are rounded off 
when Zq > 0, since the interior wave function interpolates 
smoothly between the values at z = —Zq and z = zq. A 
general ID two-body wave function ip{z) is the sum of 
even and odd parts: 'ip{z) = ipeiz) + V'o(-z)) and the zero- 
range limit of Eqs. (1) can be combined into 

^^'(0+) - V'(O-) = -(a?^)-M^(0+) + V(0-)] 
V(0+)-V(0-) = -a?^[V'(0+)+V'(0-)] . (3) 

III. EVEN AND ODD- WAVE 
PSEUDOPOTENTIALS 

Take the Hamiltonian to be 

HiD = -ih''/2^,)d^^+vtj, + v",j, (4) 

where vfjj and are even- and odd-wave pseudopoten- 

tials to be determined. 9^ is nonsingular for z 0, but 
at the origin there are singular contributions. The first 
derivative is a^VC-z) ='((>' {z^0) + [V'(0+) - 'ipiO-)]S{z). 
The second derivative then has two contributions in ad- 
dition to ^ {z 0), one because in general V' (0-I-) ^ 
V' (0— ) and the other from the derivative of the delta 
function: 

d^M^) = f{z ^ 0) + [^^'(0-^) - i^'{o-mz) 

+ [V,(0+)-V(0-)]5'(z) . (5) 

The 5{z) term is standard in the theory of zero-range 
even-wave interactions, the derivative discontinuity in 



tp^z) being chosen to cancel a zero-range even-wave in- 
teraction proportional to 6{z), but the 6 {z) term is new. 
To cancel it an odd-wave pseudopotential proportional 
to S (z) suggests itself. However, the discontimiity in 
ijj leads to unwanted products of delta functions unless 
regularizing operators are included. Define two linear 
operators i5± and d± by 

S±i^{z) = (l/2)[V>(0+) + V(0-)]5(z) 

d±7p{z) = (l/2)[V>'(0+) + V'(0-)] (6) 

where 5{z) is the usual Dirac delta function. The 
divergence-free even and odd-wave pseudopotential op- 
erators are then 

vtD=9tDS± , vIj, = glJ {z)d± . (7) 

They satisfy convenient projection properties v^jj^Jq = 
Vij^tpe = on the even and odd parts of ip, and their 
matrix elements arc (xlvinli^) = ^X*(0)[V'(0+) + V'(O-)] 
and {xKdW = -i[x'(Or[V^'(G+) + ^'(0-)]. They are 
unambiguous and connect only even to even and odd to 
odd wave functions if we stipulate that x(0) = [the 
average of x(0+) and x(0— )] if x is odd and x (0) = 
[the average of x (0-I-) and x (0—)] if X is even. In fact, 
the wave function and its derivative at 2; = refer to 
the internal wave function as modified by the potential, 
whereas z = 0+ and z = 0— refer to the wave function 
just outside the range of the potential, and the above 
values a.t z = follow from the way the internal wave 
function interpolates between the contact conditions on 
the exterior wave function (see below). Terms in S{z) 
and 5 {z) cancel from Hid if 

5m[V'(0+)+V(0-)] = (;iVM)[V''(o+)-v'(o-)] 
9iDbP'{o+) + 7p'{o-)] = (;iVM)[V'(o+)-V(o-) (8) 

and these are equivalent to the contact conditions (3) if 

9iD = -^Vmoizj and g^^ = -^^^idIiJ- 

The physical significance is clarified by starting from 
a nonsingular square well. Take the potential v{z) to 
be — Vo when —zq < z < zq and zero when \z\ > zq. 
(The odd-wave interaction i^ -^10 is negative def- 
inite in the regime of interest, where g^j-, > 0.) The 
antisymmetric solution iIj,, of the zero-energy scattering 
equation [(— ?i^/2^)a^ + v{z)]y:o{z) = inside the well 
is sin(fi;z) with k = ^/2jjVq/W. The odd-wave scat- 
tering length a^jj is defined by the second Eq. (1), 
which is satisfied in the limit zq 0-|- if scales 
with zo as k = (7r/2zo)[l + (2/7r)2(zo/a°^)]. In that 
limit the boundary conditions reduce to the second 
Eq. (3). Inside the well the kinetic and potential en- 
ergy terms are —{t^/2y)dlipo{z) = — (fi.^«;^/2y(i) sin(Kz) 
and v{z)'4'o{z) = — Vb sin(«;z). For jzj < zq, cos(kz) 
is proportional to a representation of 5{z) as zq — > 
0, since /^"^ cos{nz) f {z)dz /(O) /^°^ cos{Kz)dz = 

/(0)2k~^ sin(KZo) 2zof{0)- Then its derivative 
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— (k/2zo) sin(Kz) is a representation of 6 (z). Noting 
that kzq — > 7r/2 as zq — > we have —{h"^ /2ij.)d^ipo{z) = 

— {h?K^/2^)sm{Kz) — > {nh'^ /2^)6 (z) which agrees with 
the kinetic energy term — (fi,^/2y(i)[V'o(0+) — tpo{0—)]S (z) 
from Eq. (5) since ■0o(O+) and Vo(0— ) are to be inter- 
preted as ipo{zo) and Vo(— -^o) as zq — >■ 0+. Next con- 
sider the potential energy term inside the well as zq — > 
0+: -Vosm{Kz) - Vo(-2zo/ K)d' (z) {-nh^ /2iJi)6' {z). 
Comparing this with v°jjipo{z) from Eq. (7), using the 
expression for g°jj, noting that tp (0±) in Eq. (6) are 
to be interpreted as ip {±zo), one finds that the two ex- 
pressions for the potential energy term agree in the limit 
zo 0+. 



IV. FERMI-BOSE MAPPING 

The two-body states ipi^) considered so far are 
fermionic, i.e., the spatially even part V'e(z) contains 
an implicit spin-odd singlet spin factor, and the spa- 
tially odd part ipo{z) contains implicit spin-even triplet 
spin factors. To emphasize the combined space-spin 
fermionic antisymmetry, these will now be denoted by 
iPf{z) = i>p{z) + t/jpiz). States of combined space- 
spin bosonic symmetry can be defined by the mapping 
tpB{z) = sgn{z)tpp{z) where sgn(z) is +1 if z > and 
— 1 if 2 < 0. This maps the spatially even fermionic 
function tj^p to a spatially odd bosonic function tpg 
and the spatially odd fermionic function to a spa- 
tially even bosonic function ipg while leaving the spin 
dependence unchanged, and the corresponding scatter- 
ing lengths are also unchanged: alp, p = a\p p and 
a\p p = a°p p. Then the even- wave contact conditions 
for a^£) p follow from the odd-wave contact conditions 
for a^p p and the odd-wave contact conditions for a\p p 
follow from the even- wave contact conditions for a\p p. 
Since the kinetic energy contributions from z ^ Q also 
agree, one has a mapping from the fermionic to bosonic 
problem which preserves energy eigenvalues and dynam- 
ics. The bosonic Hamiltonian is of the same form as 
the fermionic one (4) but with mapped coupling con- 
stants glp p = h^ii^glp p and g^p p = K^/ fi'^gfp p, 
the first of which agrees with the low-energy limit of 
Eq. (25) of [9, 19]. In the limit gfpp = +oo arising 
when Vp Q—, this is the N = 2 case of the original 
mapping [10, 11] from hard sphere bosons to an ideal 
Fermi gas, now generalized to arbitrary coupling con- 
stants and spin dependence. This generalizes to arbi- 
trary N: Fermionic solutions ippizi, ui; • • • ; Zn, cjv) are 
mapped to bosonic solutions ipBizi, (Ji; • ■ ■ ; zj\j, (7j\j) via 
TpB = A{zi, - ■ ■ , ZN)ipF{zi,ai; - ■ ■■,ZN,a-N) where A = 
Yli<j<£<N ^S^izje) is the same mapping function used 
originally [10, 11] and the spin z-component arguments 
aj take on the values t and J,. The N-fermion and 
N-boson Hamiltonians are both of the form Huj = 
i<j<e<N 

generalizing (4) and (7), where the linear operators 



Sji and djf are defined on the Hilbcrt space of N- 
particle wave functions ip by Sjitp — {l/2)['ip\zj=zt+ + 
ip\z,=zt-]S{zj - Zi) and djiip = {l/2)[dz^'>p\z,=zi+ - 
dzfijj\zj=zi-]- On fermionic states ipp, gfp and g°p 
are g^p p and g°p p, whereas on the mapped bosonic 
states Vb = Atpp they are g^p p = ^"'/m^Sidf ^^'^ 

V. N-PARTICLE GROUND STATE 

Assume that both gfp p > and g^p p > 0. If g^p p 
is zero or negligible, then it follows from a theorem of 
Lieb and Mattis [22] that the fermionic ground state has 
total spin S = (assuming N even), as shown in the 
spatially uniform case by Yang [23] and with longitudi- 
nal trapping by Astrakharchik et al. [24] . If g^p p is not 
negligible then the ground state may not have 5 = 0. 
In fact, if glp p is zero or negligible then one can ap- 
ply a theorem of Eisenberg and Lieb [25] to the mapped 
spinor boson Hamiltonian, with the conclusion that the 
degenerate Bose ground state is totally spin-polarized, 
has S = N/2, and is the product of a symmetric spatial 
wave function i/»bo and a symmetric spin wave function. 
ipBo is the groimd state of the Lieb-Liniger gas [14, 15] 
which is known for all positive gljj B (hence all posi- 
tive mapped g°p p) in the absence of longitudinal trap- 
ping. The inverse mapping then yields the A'^-fcrmion 
ground state, which has a totally spacc-antisymmctric 
and spin-symmetric wave function, which is {N + l)-fold 
degenerate since Sz ranges from —N/2 to N/2. Define 
dimensionless bosonic and fermionic coupling constants 
by 7b = mgfp p/nfi^ and jp = mg°p pTi/H^ where n 
is the longitudinal particle number density. They sat- 
isfy jBjF = 4. The energy per particle e is related to 
a dimensionless function 6(7) available online [26] via 
e = (/i^/2m)n^e(7) where 7 is related to jp herein by 
7 = 7b = 4/7i?. This is plotted as a function of -yp in 
Fig. 1. For g°p ^ ^ 00 as occurs at a p-wave Feshbach 
resonance, one has a "fermionic TG gas" [13] mapping to 
a zero-energy ideal Bose gas, a fermionic analog of the 
"TG gas" of impenetrable point bosons mapping to an 
ideal Fermi gas [1, 3, 7, 10, 11, 27]. Any S = state 
has a higher energy in this case; in fact, for A'' > 2 the 
mapped Bose gas is partially space-antisymmetric, rais- 
ing its energy by the exclusion principle. 

In the presence of a uniform external magnetic field 
h the directional degeneracy is lifted and the above TV- 
particle state is the (now nondegencratc) ground state 
with field quantization direction parallel to the field, and 
the ground state energy is lowered by an amount N^Bh/2 
where HBh/2 is the magnetic moment of each spin-i 
atom. 

So far we have considered only the extremes of an 5 = 
ground state (large gfp p) or one with S = N/2 (large 

g°P p). The determination of the state of lowest energy 
for arbitrary values of these coupling constants is as yet 
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FIG. 1: Log-log plot of scaled ground state energy per particle 
e = 2mt/h^n^ for the spatially antisymmetric spinor Fermi 
gas, versus dimensionless fermionic coupling constant 7f. 

only partially solved, although we have recently obtained 
exact results for the phase diagram of ground-state total 
spin [28]. 

VI. RESPONSE TO A MAGNETIC FIELD 
GRADIENT 

Suppose that the spinor Fermi gas is longitudinally 
trapped by an optical potential Vtrap ~ |"^^iong^j i 
and that there is also a longitudinal magnetic field 
h{z) = cz with constant gradient c, adding an interac- 
tion term Vspacc-spin = -^J■BcJ2j ^jzZj to the A^'-particle 
Hamiltonian H[]j , where Sjz is the spin z-component op- 
erator for the jth particle. The space-spin interaction 
terms can be eliminated by a canonical transformation 
U~^ZjU = Zj — asjz, U~^PjU = pj, U~^SjzU = sj^ 
which leave the canonical commutation relations invari- 
ant. Noting that ?7~Vspace-spinC^ = Vspe,cc-spin + jN^lBCa 

and U^'^VtiapU = Y,j 5"^-^iong(^j ^"^i^)^' fi^^ds that 
the space-spin coupling terms cancel from U~^H^j-,U 
with the choice a = —fi-Qc/{mLo'^^^^, leading to a 

transformed Hamiltonian U~^HijjU = Hijj{c = 0) — 



iV(/UBc)^/(8ma;iQjjg) where H^jy{c = 0) does not include 

V^space-spin- The grouud state ^boI-^i, si; ■■■ ,zn, sjv) of 
(j-'^H^ijU is the same as that of 7?^^(c = 0). The cor- 
responding single-particle density no(z) is centered on 
z = 0. It is not known analytically in the presence of lon- 
gitudinal trapping, but accurate numerical results have 
been calculated by a local density method [16]. In the 
presence of Kspacc-spin the single-particle density n{z) is 
{<pBo\U^^h{z)U\(pBo) where n{z) = J2j - zj). The 
state \(l)Bo) is a simultaneous eigenstate of the longitu- 
dinal spin operator Sz = X]j ^^jz-, which has eigenvalues 
Sz = -^N,-^N + 1,--- -1,^N. The ground state 

of Hiu, which now includes Vspaco-spin, is U\(f>Bo}, and the 
(7V-|- l)-fold degeneracy is not lifted by the magnetic field 
gradient so long as the magnetic field vanishes at z = 0. 
One has U~^6{z - Zj)U = 5{z - Zj - J^Jl" sjz), whose 

expectation value is N~^no{z =F ^^o) when the eigenvalue 
of Sjz is ±|, where Zq = ^ . If the eigenvalue of Sz 

long 

is Sz then wN of the Sj^ have eigenvalue ^ and (1 — w)N 
have eigenvalue — ^, where the fraction w, which satis- 
fies < w; < 1, is w; = N-'^Sz + It follows that the 
single-particle density n{z) in a ground state U\4>Ba) with 
longitudinal spin Sz is a weighted average of the extremal 
densities: n{z) = wnoiz — zq) + (1 — w)n(){z + Zq). The 
ground state wave function of the corresponding spinor 
Fermi gas differs by a factor A, the previously-given map- 
ping function. It has the same longitudinal spin eigen- 
values and same degeneracy, and since = 1 these 
fermionic ground states have the same density profiles 
as the bosonic ones. 
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